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We have obtained the first estimates of the masses of the components of the Her X-l/HZ 

Her X-ray binary system taking into account non-LTE effects in the formation of the H 7 

absorption line: m x = 1.8M and m v = 2.5M . These mass estimates were made in a Roche 

model based on the observed radial- velocity curve of the optical star, HZ Her. The masses for 

ly-) . the X-ray pulsar and optical star obtained for an LTE model lie are m x = O.85±O.15M and 

m v = 1.87 ± O.13M . These mass estimates for the components of Her X-l/HZ Her derived 

from the radial-velocity curve should be considered tentative. Further mass estimates from 

high-precision observations of the orbital variability of the absorption profiles in a non-LTE 

\ model for the atmosphere of the optical component should be made. 
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INTRODUCTION 



The X-ray source Her X-l was detected by the Uhuru satellite in 1972 pQ and identified 
with the variable star HZ Her [21 [3J, IH E] . The Her X-l/HZ Her binary system consists 
of the A7 optical component (6], which fills its Roche lobe, and an X-ray pulsar with a 
period of 1.24 s. The orbital period of the binary is Porb = 1.7d [7]. The system also 
displays 35-day X-ray variability due to the precession of the warped accretion disk [HE]- 
The optical variability of the binary is mainly due to X-ray heating of the optical star 

The Her X-l/HZ Her binary has been thoroughly studied: the high-precision half- 
amplitude of the radial-velocity curve for the X-ray pulsar K x = 169.049 km/s [9] has 
been derived from X-ray timing, and the orbital inclination i = 81° — 88° P, [H)] from the 
duration of the X-ray eclipses. However, no unambiguous estimates for the mass of the 
X-ray pulsar have been obtained, due to the very intense X-ray heating, k x ~ 150. 

The value of a neutron star in the Her X-l/HZ binary is not determined unam- 
biguously. According to Middleditch and Nelson [11], the mass of the X-ray pulsar is 
m x ~ 1.30 ± 0.14M Q . Koo and Kron [12] suggest that the mass of the compact object 
in Her X- 1 is close to m x ~ 1.5M . The mass estimated by Hutchings et al. [32] is 
m x = 0.93 ± 0.07M©, while Reynolds et al. [10] conclude this mass lies in the inter- 
val m x = 1.5 ± 0.3M Q . These estimates were made using point-mass models, except 
for the estimate of Koo and Kron [12], which was made for a Roche model assuming a 
spherical shape for the optical component and roughly taking into account heating of this 
component. 

The difficulty in accurately estimating the pulsar mass is essentially due to the intense 
X-ray heating of the optical component. The X-ray luminosity of the compact object 
is L x = 3.39 • 10 37 erg/s [6], causing the spectral type of the optical star to vary from 
A7 (<f) = 0.0) to B3-B6 (<fi = 0.5) over half the orbital cycle [61 [H]. The temperature 
distribution over the surface of the optical star is also complicated by shielding of the 
X-ray radiation by the warped, precessing accretion disk. Together, these factors hinder 
efforts to measure a high-precision, purely orbital radial velocity curve, and thereby derive 
an unambiguous estimate for the mass of the X-ray pulsar. 

Note that, in a point-mass model, only the halfamplitude of the radial- velocity curve, 
K v , is taken into account, not its shape. However, it is known that the observed radial- 
velocity curves of stars in binary systems may differ from the radial-velocity curves for 
the centers of mass, due to tidalrotational distortion of the stellar surfaces and to heating 
by the radiation from their companions (see the studies of Wilson and Sofia [15], Milgrom 
[T6] , and references therein). It is particularly important to take these effects into account 
when determining themasses of relativistic objects in X-ray binaries. For example, Milgrom 
[T7] [To] showed that, even in the case of relatively low X-ray heating of the optical star 
(as in Cyg X-l and SMC X-l), the radial velocities derived from absorption lines can 
vary substantially. Using profile calculations for various absorption lines in the spectra 
of HZ Her and other stars [T7J [16], he noted the important role of the "soft" X-ray 
radiation (E < 1 keV) in the spectrum of the relativistic object, which is responsible for 
the emission components of the lines. Milgrom and Salpeter [TS] developed a procedure 
for calculating the radiation emerging from the atmosphere of a star irradiated by X-ray 
radiation, intended for the calculation of line profiles. Unfortunately, in the 1970s, such 
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calculations required enormous amounts of computer time, and could not be applied to 
the massive calculations required to solve this problem. 

Antokhina et al. [3Ql El] suggested a procedure for calculating line profiles and radial- 
velocity curves of optical stars in X-ray binaries taking into account the incident X- 
ray flux from the relativistic object. Using this technique, we carried out a series of 
studies analyzing spectroscopic observations of close X-ray binaries [19j |20j [211 EZ], and 
demonstrated the need to take into account the shape of the radial-velocity curve when 
determining the component masses, particularly in the presence of intense X-ray heating. 

Here, we analyze the observational data of pi)] in a Roche model exactly taking into 
account the effect of X-ray heating of the optical star [281 [291 E] , with the aim of reliably 
estimating the mass of the compact object in Her X-l/HZ Her. Running ahead for a 
moment, we note that no emission components are visible in the absorption lines in the 
spectrum of the optical star HZ Her, despite the extremely intense X-ray heating. Since 
this may be related to deviations from LTE, we have analyzed the radial-velocity curve 
using both LTE and non-LTE models. 

OBSERVATIONS 

Let us briefly describe the spectral data adopted from [TO]. The observations were 
carried out on June 10-16, 1995 with the 2.5-m Isaac Newton Telescope of the Royal 
Greenwich Observatory at La Palma (Canary Islands). In total, 59 spectra at wave- 
lengths of 4080-4940 AA and with exposures from 1200-1800 secunds were obtained with 
the Intermediate Dispersion Spectrograph equipped with a 235-mm camera, a diffraction 
grating with a dispersion of 1200 lines/mm, and a 1024x1024 pix Tektronix CCD array. 
The inverse linear dispersion of the spectrograms was 0.84 A/pixel. The signal-to-noise 
ratios in the spectra obtained were S/N = 50 — 100 The wavelength scale was calibrated 
with a Cu-Ar standard; the calibration error was < 0.05 A. 

The radial velocity of HZ Her was determined via a cross-correlation analysis us- 
ing the radial velocities of IAU standard F stars observed on the same nights as HZ 
Her. The epoch of the maximum positive radial velocity of the optical component 
(JD=2448800.537) was used as the zero phase in [10]. We adopted the orbital period 
P orb = 1^.700167412 0. 

Here, we took the zero orbital phase to correspond to the middle point of the eclipse of 
the X-ray source by the optical component. The data analysis of [TO] shows that, in this 
case, the maximum positive radial velocity of the optical star falls at phase <p — 0.22, and 
the systematic radial velocity of the center of mass is 7 = —69 km/s. Figure [Q presents 
the observed radial-velocity curve. 

The spectrograms of [10] were obtained for the orbital-phase interval of the 35-day 
precessional cycle of the accretion disk 35 = 0.14 — 0.33 (the zero phase was ascribed to 
the moment when the X-ray source "switches on" for an observer on Earth). Thus, at 
the epoch of the spectroscopical observations [10J , the accretion disk did not shield the 
optical star from the X-ray flux, and our fitting of the observed radial- velocity curve did 
not take into account this shielding. 
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MODEL OF THE BINARY AND PROCEDURE FOR CALCULATING THE 
RADIAL- VELOCITY CURVE 



We synthesized theoretical absorption-line profiles and radial-velocity curves for the 
optical star using the two procedures described in detail in [281 123 EI], which we call 
Procedure I [28j[29] and Procedure II [30j[31]. These procedures are similar, and differ only 
in the means used to calculate the radiation flux from local areas of the star. Procedure 
I (the more rapid algorithm) uses the Balmer line profiles tabulated by Kurucz [32] ; 
in Procedure II (the more accurate algorithm), the emergent flux is calculated using a 
model atmosphere in the presence of incident external flux. We will briefly describe the 
calculation technique. 

In a close binary, the optical component has a tidally distorted shape and an inho- 
mogeneous temperature distribution over its surface, due to the effects of gravitational 
darkening and surface heating by the X-ray radiation from the relativistic object. There- 
fore, to model the radial-velocity curve, the X-ray binary was represented by the optical 
star in the Roche model and an X-ray point source with a finite mass. The tidally de- 
formed stellar surface was divided into ~ 2600 elements, for each of which the emergent 
radiation was calculated. The calculations of the flux from a surface element included 
gravitational darkening, the heating of the stellar surface by radiation from its compan- 
ion (the reflection effect), and limb darkening. In Procedure I [28, 29J, the heating of 
the star's atmosphere by X-ray radiation from the companion was calculated by adding 
the emergent and incident flux without considering radiative transfer in the stellar atmo- 
sphere. 

For each visible surface element with temperature Ti oc and local gravitational acceler- 
ation gi oc , the absorption-line profiles and their equivalent widths were interpolated from 
the Kurucz tables for Balmer lines [32] . Summing the local profiles over the visible surface 
of the star, taking into account the Doppler effect, and normalizing the profiles to the con- 
tinuum level for each surface element, we calculated the integrated profile from the star at 
a given orbital phase (see [2S1 US\ for details). The calculated integrated absorption-line 
profile was used to determine the radial velocity of the star. The radial velocity at a 
given orbital phase was calculated for the average wavelength at the levels of the residual 
intensities at 1/3, 1/2, and 2/3 of the integrated absorption profile. 

In the more modern Procedure II, for each surface element, the parameter k l ° c , equal 
to the ratio of the incident X-ray flux and emergent radiation flux, is calculated separately 
from the local temperature T} oc and local gravitational acceleration gi oc , without allowing 
for the external irradiation of the atmosphere. The model atmosphere is calculated for 
these parameter values at a given point of the surface, assuming LTE, by solving the equa- 
tions of radiative transfer in the line in the presence of incident external X-ray radiation 
(see [3UJ EI] for details). Thus, the intensity of the emergent radiation in the line and 
continuum is calculated for each local surface element. Further, following a procedure 
identical to Procedure I, the integrated absorption-line profile at the given orbital phase 
is calculated, which is then used to determine the radial velocity of the optical star. 

In summary, the main difference between Procedure II [33 EI] and Procedure I [2"8"| |2"5] 
is that, in the former the local profile of a surface element is found by constructing an LTE 
model of the atmosphere and calculating the intensities of emergent radiation in the line 
and in continuum (taking into account reprocessing of the external X-ray radiation). In 
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contrast, Procedure I uses calculated and tabulated Kurucz Balmer absorption-line profiles 
for various effective temperatures T e ff and gravitational accelerations g. Note that, in 
addition to the simplified treatment of the reflection effect, this method for calculating 
the hydrogen absorption profiles is not entirely correct, since the theoretical line profiles 
in Kurucz's tables [32] are given in relative fluxes rather than intensities. However, since 
we use the theoretical line profiles to determine radial velocities and not for comparison 
with the observed profiles, we consider this approximation to be reasonable [21 J . 

Substantial deviations from LTE can occur in the atmospheres of hot stars, particularly 
in the presence of external irradiation. Therefore, we modeled the Balmer line profiles 
with Procedure II both assuming LTE and taking into account non-LTE effects in the 
hydrogen atom. The equilibrium non-LTE populations of the HI levels were calculated 
using a technique similar to that developed by Ivanova et al. [22]. In the calculations, 
a 23-level model of the HI atom was applied, taking into account all allowed collisional 
and radiative mechanisms for the atom's redistribution over the states, including those 
generated by the external radiation. In total, the HI model includes 153 allowed bound- 
bound and 22 bound-free radiative transitions, 55 of which were linearized. The non-LTE 
populations of the levels were obtained for stellar-atmosphere models for the local surface 
elements using the NONLTE3 [23] package, which implements the total linearization 
technique of Auer and Heasley [24J, modified by Ivanova et al. [22] to take into account 
the external radiation. The calculations included all continuum sources of opacity at 
optical and X-ray wavelengths and about 570 000 lines from the lists of Kurucz [25] . 
The subsequent modeling of HI line profiles was carried out following the above standard 
technique, but using the non-LTE populations. 

When the local profiles in FigfS] obtained for the LTE and non-LTE models are com- 
pared, it is obvious that the H 7 line does not display an emission component in the 
non-LTE case. This is due to the strong (by a factor of four to six) overpopulation of the 
second level (n = 2) of the hydrogen atom (relative to the above levels) in the chromo- 
sphere. The source function S u ~ j-B v (T e ) is proportional to the ratio of the populations 
of the lower and upper levels of the transition. As a result, in spite of the substantial 
increase of B u (T e ) in the hot chromosphere, S u exceeds the source function in the region 
of formation of the continuum; i.e., the line profile appears in absorption. 

As follows from an analysis of physical processes in the HI atom, an overpopulation 
of its low-excitated states should be characteristic of hot objects with Xray irradiation. 
According to [26], hydrogen is responsible for less than 20% of the total absorption at 
X-ray wavelengths for T e = 7000-fT, and for a factor of 10 3 less for T e = 20000-fT. 

At the same time, since the temperature in the chromosphere exceeds that in the region 
of continuum formation, the source function in the chromosphere dominates over the 
mean intensity of the radiation at all optical frequencies, J v <C S u , where J v is the mean 
intensity of the radiation and S v the total source function, both at a given frequency at a 
given depth. Thus, radiative recombination to all HI levels and subsequent spontaneous 
transitions, which are proportional to S v , dominate over ionization processes due to both 
the external radiation and intrinsic radiation of the atmosphere. As a result, all HI levels 
turn out to be overpopulated relative to the LTE approximation, with the overpopulation 
increasing for lower states. 

These non-LTE effects should result in a strengthening of the absorption and weak- 
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ening of the emission components of hydrogen lines in the spectrum of an atmosphere 
with substantial X-ray irradiation. We wish to emphasize that, in the spectrum of the 
optical star HZ Her, in spite of the huge amount of X-ray heating (k x ~ 150), no emission 
components are observed in the profiles of absorption lines, possibly due to deviations 
from LTE. 



ANALYSIS OF THE OBSERVED RADIAL- VELOCITY CURVE 

We fitted the observed radial- velocity curve of the optical star within the Roche model. 
Table [1] presents the model parameters for the X-ray binary Her X-l/HZ Her. 



Table 1: Parameters used in the modeling of the radial- velocity curves of the optical 
component in the Roche model. 




1.700167412 
0.0 



Orbital period 
Eccentricity 
Orbital inclination 

Degree of Roche lobe filling by the optical component 
Asynchronicity factor for the rotation of the optical compo- 
nent 

Effective temperature of the non-perturbed optical compo- 
nent 

Gravitational-darkening coefficient 

Ratio of the X-ray luminosity of the relativistic component 
and the bolo metric luminosity of the optical component, 
L x / L v 

Coefficient for reprocessing of the X-ray radiation (Procedure 
I) 

Limb-darkening coefficient 



As is noted above, Her X-l/HZ Her displays strong X-ray heating of the atmosphere of 
the optical star. The X-ray luminosity of the compact object is L x ~ 3.39 ■ 10 37 erg/s [6], 
whereas the luminosity of the optical star is L v ~ 2.28 • 10 35 erg/s (the optical luminosity 
was obtained in the Roche model assuming the mean effective temperature T e ff = 8100K 

ED- 

In our modeling of the heating of the atmosphere of the optical component, we used 
the X-ray spectrum of Her X-l obtained with BeppoSAX in 1998 [33]. In spite of the 
strong interstellar absorption, soft X-ray radiation from Her X-l at 0.1-1 keV is reliably 
detected, which is essential to study the conditions for the formation of the optical spec- 
trum of HZ Her [17]. For example, Fig ]3] presents two model LTE spectra of the H 7 
hydrogen absorption line of the optical component calculated with and without taking 
into account the soft component (0.1-1 keV) of the irradiating X-ray spectrum. When 
this soft component is taken into account, the emission component in the H 7 absorption 
line is appreciably stronger. 

The X-ray spectrum at 0.1-10 keV, corrected for the interstellar absorption, was re- 
constructed from the data of [33] using the XSPEC code. The two-component model of 
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the spectrum was taken from [33]: blackbody radiation with the temperature 0.093 keV 
(which dominates up to ~ 1 keV), and power-law-like radiation with spectral index 0.74 
(which dominates at energies exceeding 1 keV). We did not include the components from 
[33] that describe the X-ray iron emission lines due to their small contribution to the 
total X-ray flux. Whenmo deling the optical spectrum of HZ Her, we assumed the X-ray 
radiation of the relativistic object to be isotropic. 

Since the exact mass of the optical star in the Her X-l/HZ Her system is not known, 
we adopted the masses of both components parameters to be fit. We obtained 

the fit via an exhaustive search over the parameters, making it possible to study the 
surface of residuals in terms of various parameters in detail. 

The half-amplitude of the radial-velocity curve of the X-ray pulsar, K x = 169.049 
km/s [9], is known from X-ray timing, and consequently the mass function of the pulsar 
is f x (m) = O.85M . Therefore, in the course of the exhaustive search over the component 
masses, we particularly made sure that the mass function of the X-ray pulsar retained the 
value f x {m) = O.85M . This was done as follows. 

The mass function 

m^sin 3 i 

fx{m) - 



(m x + m v ) 2 ' 

for a given orbital inclination i specifies an unambiguous relation between the masses of 
the optical star m v and the X-ray pulsar m x . 

A value of m x was calculated from ([T]) for each value of m v and a fixed value of i. 
Further, the theoretical radial-velocity curve for the optical component was determined 
based on the resulting component masses m x and m v , with the numerical parameters of 
the Roche model taken from Table [H The adequacy of the model to the observational 
data was tested using the statistical \ 2 criterion, selecting a significance level of 5% (see 
for details). 



Results Obtained with Procedure I 



Since the radial velocities of [10] were derived from spectra obtained with the response 
function FWHM=1.7A, we calculated the theoretical radial- velocity curves from the pro- 
file of the H 7 line convolved with an instrumental profile with FWHM=1.7A. Figure H] 
presents the theoretical integrated profiles of the H 7 absorption line. 



Table 2: Mass of the relativistic and optical components of Her X-l obtained in the Roche 
model with Procedure I. 



Orbital 
inclination 


m x (M Q ) 


m v (M Q ) 


i = 80° 
i = 88° 


1.84 
1.78 


2.6 
2.5 



Table [2] presents the masses of the X-ray pulsar m x and optical star m v corresponding 
to the minimum \ 2 residual for the theoretical and observed radial-velocity curves. In 
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Table [2] the component masses are given without errors, since models of the binary were 
rejected according to the selected significance level a = 5%. For the quantile A (a = 5%) = 
77.93, the minimum residual for the orbital inclinations i = 80°, i = 88° was xLin — 170. 
Figure [5] presents the theoretical radial- velocity curve and Figj6k the behavior of the \ 2 
residual for i = 88°. 

Results Obtained with Procedure II 

It was noted above, we used the observational data of \±0\ ■ In those spectra, the profile 
of the H 7 absorption line at orbital phase 0.47 displays a distortion [[10], Fig. 1]. At the 
same orbital phase, the theoretical integrated profile of the H 7 absorption convolved with a 
response function with FWHM=1.7 A observed profile (Fig. [7]). We were able to reproduce 
a similar distortion in the integrated theoretical profile of the H 7 hydrogen line at phase 
0.47 only with a response-function width of FWHM=5 A (Fig. [7j). Therefore, we carried 
out our analysis of the observed radial-velocity curve for the H 7 hydrogen line with the 
response-function width FWHM=5 A. 

For the effective temperature of the stellar surface ~ 8000 K, the hydrogen line displays 
a full width at half maximum of the residual intensity of several tens of A. In our case, 
the width of the model integrated profile of the H 7 line at this level was ~ 25 A (Fig. 
IS}. Therefore, we faced the problem of determining the radial velocity from such a broad 
H 7 absorption profile. For Procedure I, the half-width of the model profile was smaller 
~ 10 A (Fig. IS]), and the profile itself was more symmetrical. Therefore, to determine the 
position of the center of gravity of the line, it was sufficient to take three cross sections 
at heights of 1/3, 1/2, and 2/3 of the residual intensity. If the integrated H 7 profile is 
calculated with Procedure II, three cross sections are not sufficient, and we selected six 
and eight levels of the cross section of the residual intensity of the model integrated profile 
to determine the center of gravity of the H 7 line. 

Table 3: Masses of the relativistic and optical components obtained from the radial ve- 
locity determined from six cross sections of the integrated H 7 . 



Orbital 
inclination 


m x (M Q ) 


m v (M Q ) 


i = 80° 
% = 88° 


0.85±0.15 
0.81±0.13 


1.87±0.13 
1.80±0.11 



In the former case, the radial velocity was determined from the wavelength averaged 
over cross sections at six levels of the lower half of the absorption profile; the resulting error 
in the radial velocities were ~ 1 km/s. When eight cross-section levels were used, two cross 
sections in the upper half of the model profile of the H absorption line were introduced 
in addition to the six for the previous case; the resulting error in the theoretical radial 
velocities were ~ 2.5 km/s. Here, we take the radial velocity at orbital phase <fi = 0.0 
(when the optical component eclipses the X-ray source) as the error in the theoretical 
radial velocity. The error in the radial velocity for the calculation using six cross sections 
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is appreciably smaller due to the narrowness of the lower part of the H 7 profile compared 
to the wings (Fig. EJ). 

The mass of the X-ray pulsar m x and optical star m v corresponding to the minimum 
X 2 residual for the theoretical and observed radial-velocity curves obtained using six and 
eight cross sections of the integrated H 7 line profile are presented in Tables [3] and HJ 
respectively. Figure H 7 presents the radial- velocity curve obtained using six cross sections 
of the integrated |5] profile. 

It is obvious that the results depend on the method used to determine the center of 
gravity of the integrated H 7 absorption profile, although, within the errors, the mass of 
the X-ray pulsar is consistent with that of the optical star. 

Table 4: Masses of the relativistic and optical components obtained from the radial ve- 
locity determined from eight cross sections of the integrated H 7 profile. 



Orbital 
inclination 


m x (M Q ) 


m v {M Q ) 


i = 80° 
i = 88° 


0.70 
0.70 


1.75 
1.70 



The model of the binary obtained using six cross sections is consistent with the ob- 
servational data at the 7 = 95% confidence level, with the minimum residual Xmin — ^ 
(Figj6b). When eight cross sections are used, the binary model is rejected at the signifi- 
cance level a = 5%, with the minimum residual Xmin — 90 (f° r this reason, the component 
masses in Table H] are presented without errors). 

Integrated Profile and the Radial-velocity Curve Including non-LTE Effects 

As is noted above, in spite of the substantial X-ray heating (k x = 150), the absorption 
lines in the optical spectrum of the HZ Her system display no visible emission components. 
This is apparently due to non- LTE effects (see above). 

We also calculated the integrated H 7 absorption profiles and corresponding radial- 
velocity curve in a non-LTE model (Fig. [5]). These calculations were performed using 
the NONLTE3 code [23] for the model binary parameters m x = 0.81M Q , m v = 1.8OM , 
and % = 88°. (see the remaining parameters in Tabled]), taking the width of the response 
function to be FWHM=1.7A. Regrettably, calculations with the NONLTE3 procedure 
are time consuming, which currently makes searches for the parameters of close binaries 
challenging, and we restrict our consideration here to qualitative conclusions. 

Figures [2] and |H] present the local and integrated profiles of the H 7 absorption line 
calculated in the non-LTE approximation. Figure E] presents the radial- velocity curve 
derived from these profiles. 

We can see from FigJH]that, qualitatively, the integrated H 7 absorption profile obtained 
in the non-LTE model is close to the integrated profile calculated using Procedure I 
- in both profiles, no emission component is seen. At the same time, the integrated 
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H 7 absorption profile obtained assuming LTE with Procedure II displays an emission 
component (see Figsj3] and E]). 

We can also see in Figj5] that the shape of the radial- velocity curve obtained with 
the non- LTE model is qualitatively closest to that calculated with Procedure I. In this 
connection, the results of fitting the observed radial-velocity curve for HZ Her/Her X-l 
made with Procedure I may bemore reliable: m x = 1.78M Q , m v = 2.5M Q for the orbital 
inclination i = 88°. We stress again, however, that the model of the binary system with 
these parameters was rejected at the a = 5% significance level adopted in the study. 

DISCUSSION 

Let us consider the advantages and disadvantages of the procedures used to calculate 
the radial-velocity curves in our study. 

When the integrated profile of the optical component is calculated with Procedure I, 
we use integrated profiles of absorption lines of single stars tabulated by Kurucz [32j for 
the local profiles. These absorption profiles do not contain emission components, and the 
local profiles used in Procedure I are qualitatively consistent with the local profiles of the 
H 7 absorption lines in the non-LTE approximation. However, the model integrated H 7 
profile obtained with Procedure I does not reproduce the distinctive distortion of the H 7 
absorption line at orbital phase <ft — 0.47 (|10|. Fig. 1), and the best-fit theoretical radial- 
velocity curve seems to diverge from the observed radial-velocity curve both qualitatively 
and quantitatively (xLin — 170 - the model is rejected at the a = 5% significance level; 
Figj5]). In this light, we cannot take our estimates for the Her X-l/HZ Her component 
masses obtained with Procedure I to be final. 

Let us consider the results of our fitting of the observational data and the specific 
features of calculating the model absorption profiles with Procedure II. 

The local profiles assuming LTE with exact account of the effect of heating contain an 
emission component (which is not seen in the observed optical spectra of HZ Her), whereas 
the local H 7 absorption profiles calculated in the non-LTE approximation display no 
emission component (FigJ2J). Nonetheless, the theoretical integrated H 7 absorption profile 
convolved with an artificially broadened response function with FWHM=5 A (exceeding 
the real width of the response function by a factor of three) qualitatively reproduce the 
"distortion" of the observed absorption profile shape (cf. Fig. 1 in [10] and FigJS]). The 
full width of the obtained integrated H 7 absorption profile at half the residual intensity 
was ~ 25A, hindering determination of the center of gravity of the absorption line. Six or 
eight cross sections were used to determine the center of gravity of the model integrated 
H 7 absorption profile. 

The component masses obtained with six and eight cross sections of the model inte- 
grated H 7 profile are consistent within the errors (Tables [3] and H]). However, both the 
quantitative estimates of the mass of the binary components obtained with Procedure 
II and the reliability of the binary model depend on the number of cross sections of the 
integrated profiles used to determine the center of gravity of the H 7 absorption profile. 
For example, when eight cross sections of the model integrated H 7 profile were used, the 
model was rejected at selected significance level of a = 5%, and the masses for the bi- 
nary components differed by ~ O.1M from those obtained when six cross sections of the 
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integrated profile were used. 

Thus, it does not seem reasonable to give preference to any one of the various results, 
and we consider the results presented in Tables [2] and [3] equally valid. As we have shown 
here, reliable estimation of the component masses for the Her X-l/HZ Her binary system 
is possible only based on an analysis of the radial-velocity curve jointly with the orbital 
variability of the absorption profiles, with correct account for non-LTE effects. 

CONCLUSION 

Our study has mainly focused on the methodical side of the problem considered. We 
have shown that determining the radial velocities in close binaries with intense X-ray 
heating of the optical star requires taking into account non-LTE effects in the formation 
of absorption lines in the spectrum of the optical star. For Her X-l/HZ Her-type binary 
systems, not only the radial-velocity curve (which is specified by the line profiles only 
ambiguously) must be analyzed, but also line profiles at different phases of the orbital 
period. 

Our calculations have shown that, due to the very intense X-ray heating in Her X-l/HZ 
Her, the H 7 absorption line of the optical star is appreciably asymmetrical, making it diffi- 
cult to determine unambiguously the center of gravity of its profiles and, as a consequence, 
the radial velocities and component masses. For example, when the theoretical radial ve- 
locity is determined using different numbers of cross sections of the absorption-line profile, 
the mass of the X-ray pulsar m x varies by ~ 15% (Tables [3] and H]). We stress that the 
methods for estimating the Her X-l/HZ Her component masses presented here are based 
on the most current and physically justified techniques, in which local absorption-line 
profiles are calculated taking into account external X-ray heating. Nonetheless, even this 
approach cannot yield unambiguous mass estimates for the X-ray pulsar in Her X-l/HZ 
Her, due to the enormous rate of X-ray heating in this unique X-ray binary. Analyzing 
the radial-velocity curve in two ways, we have obtained two masses for the X-ray pulsar: 
m x = 0.85 ± O.15M and m x = 1.8M . Thus, the mass of the X-ray pulsar in Her X-l is 
currently determined only to within a factor of ~ 2, due to uncertainty in realistic models 
for the formation of line profiles in the spectrum of the heated optical star. 

Our results demonstrate the need to use high- resolution spectrograms (A/ A A ~ 50000) 
and to take into account non-LTE effects when estimating the mass of the compact object. 
Correct and reliable estimates for the component masses can be obtained only by analyzing 
also the orbital variability of the absorption-line profiles in the spectrum of the optical 
component in a non-LTE approximation. 

We thank I.I. Antokhin for useful discussions. 
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Figure 1: Observed and theoretical radial- velocity curves of the optical component of 
the X-ray binary Her X-l/HZ Her. The points indicate the radial velocities of HZ Her 
taken from (TO] . The solid curve is the theoretical radial- velocity curve for a Roche model 
calculated using Procedure II, with the mass of the compact object m x = O.81M , the 
mass of the optical star m v = 1.80M Q , and the orbital inclination % = 88°. (for the 
remaining parameters of the binary, see Table [T|) . 
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Figure 2: Model local profiles for the H 7 absorption line. The solid curve corresponds to 
the calculations assuming LTE (Procedure II) for m x = O.9M , m v = 2.0M Q , % = 80°, 
k x = 150, and the dashed curve to the non-LTE calculations for the same parameters of 
the binary. The temperature and gravitational acceleration for the local surface element 
are T loc = 7825K and log gi oc = 3.336. 
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Figure 3: Effect of taking into account the soft component of the incident X-ray flux on 
the shape of the H 7 line. Integrated profiles calculated with Procedure II for m x = 1M Q , 
m v = 2M Q , i = 80°, k x = 150 are shown. The X-ray irradiation was modeled at 0.1-1 
keV(dashed curve), 0.1-17 keV (dash. dotted curve), and 1-17 keV (solid curve). 
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Figure 4: Integrated profiles of the H 7 absorption line calculated in the Roche model with 
Procedure I, for the width of the response function FWHM=1.7 Aand with m x = 1.78M©, 
m v = 2.5M , and i = 88°. The solid line corresponds to the orbital phase = 0.0, and 
the dashed line to = 0.30. 
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Figure 5: Observed radial- velocity curve for the optical component of the X-ray binary 
Her X-l taken from [10] (points) and theoretical radial- velocity curves calculated with 
Procedure I for m x = 1.78M , m v = 2.5M , and i = 88° (solid curve); and with Procedure 
II for m x = 0.81M Q , m v = 1.80M Q , and i = 88° from six cross sections of the integrated 
profile of the H 7 line (dashed curve) and using the NONLTE3 [23] package in the non- 
LTE model for m x = O.81M , m v = 1.8OM , i = 88°. (dash-dotted curve). 
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Figure 6: Discrepancies between the observed and theoretical radial velocity curves of HZ 
Her. The theoretical curves were obtained with (a) Procedure I and (b) Procedure II, for 
the orbital inclination i = 80°. The horizontal line corresponds to the critical level of the 
X 2 criterion A 59 = 77.93 for the significance level a = 5%. 



18 




Figure 7: Model integrated profile of the hydrogen absorption line at orbital phase 
<p = 0.47, calculated assuming LTE with Procedure II for the response-function widths 
FWHM=1.7A,(solid curve) and FWHM=5 A, (dashed curve) for m x = O.81M , m v = 
1.8OM , and i = 88°. 
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Figure 8: Model integrated profile of the H 7 hydrogen absorption line at orbital phase 
= 0.30, obtained with Procedure II for the response-function width FWHM=5A, and 
m x = 0.81M Q , m v = 1.8OM , and % = 88° (dashed curve); with Procedure I for the 
response-function width FWHM=1.7A, and m x = 1.78M , m v = 2.5M , and i = 88° 
(solid curve); and with the NONLTE3 package in the non-LTE model for the response- 
function width FWHM=1.7A, and m x = O.81M , m v = 1.8OM , and % = 88° (dash- 
dotted curve). 
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